Helium-cluster decay widths of molecular states in beryllium and carbon isotopes 
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The a particle and 6 He emissions from possible molecular states in beryllium and carbon isotopes 
have been studied using a mean-held-type cluster potential. Calculations can reproduce well the a- 
decay widths of excited states in 8 Be, 12 C and 20 Ne. For the nucleus 10 Be, we discussed the a-decay 
widths with different shapes or decay modes, in order to understand the very different decay widths 
of two excited states. The widths of 6 He decay from 12 Be and a decays from 13 - 14 C are predicted, 
which could be useful for future experiments. 
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The experimental search for molecular-type structures 
in light nuclei is of increasing interest. The Ikcda diagram 
as a guide line reveals that cluster structures are expected 
to appear near decay thresholds [l[ ■ Evidences for dimers 
in beryllium isotopes and polymers in carbon isotopes 
have been summarized by von Oertzen [2J. It is convinc- 
ing that the candidate cluster bands should have very 
large moments of inertia, corresponding to molecular- 
like shapes. Cluster states above thresholds should be 
understandable though it is difficult to perform the pre- 
cise measurements of cluster decay widths. The model 
calculations of decay widths are generally dependent on 
the deformations and angular momenta of nuclear states. 
Hence, the predictions of widths can provide useful struc- 
ture information about tentative cluster states. 

In 8 Be, the rotational cluster band and its decay prop- 
erty have been well established based on an a- a struc- 
ture Q. In heavier beryllium isotopes, two a particles 
and additional valence neutrons can give rise to covalent 
molecular binding. Recent experiments have observed 
the a-r- 6 He decay in 10 Be d Qand the 6 He+ 6 He decay 
in 12 Be [E 0], which would be the indication for molec- 
ular states. Different molecular-type structures in 10 Be 
and 12 Be have been suggested by the antisymmetrized 
molecular dynamics (AMD) calculations [8, 9]. It is also 
intriguing that carbon isotopes with three a clusters can 
have two different configurations: triangular and linear. 
For 13 C and 14 C, cluster bands based on the a+ 9 Be and 
ct+ 10 Be systems have been proposed 0, El E3, respec- 
tively. The experimental study of a decays from excited 
states of 13 > 14 C is attracting interest [l2l.[l3l]. 

Many theoretical works based on the WKB approach 
have shown the successful calculations of decay life-times, 
e.g. 0, GJ, EE E3, El, El H3- Combining the Bohr- 
Sommerfeld quantization, Buck et al. has achieved the 
calculations of spectra of cluster states using macroscopic 
cluster potentials [2l|. Recently, we suggested a mean- 



field-type cluster potential for various charged-cluster de- 
cays from the ground states (g.s.) of even-even heavy 
nuclei (22J. In the present work, we extended our calcu- 
lations by including the deformations and angular mo- 
menta of cluster states. We focus the decay properties 
of excited states in light nuclei, particularly in beryllium 
and carbon isotopes, which are attracting great interest 
in experiments. 

The cluster potential in the quantum tunnelling ap- 
proach can be written as (e.g., (l4|). 

V(r)=V N (r) + Vc(r) + ^(L + l/2) 2 , (1) 

that contains the nuclear potential VAr(r), the Coulomb 
potential Vc{r) and the Langer modified centrifugal po- 
tential [Til- We suggested the nuclear potential Vzv(r) 



V N {r) = X[Z c v p (r) + N c v n (r)}, 



(2) 



where A is the folding factor; N c and Z c are the neu- 
tron and proton numbers of the cluster, respectively; 
v n (r) and v p (r) are the single-neutron and -proton poten- 
tials(excluding the Coulomb potential) respectively, ob- 
tained from Skyrme-Hartree-Fock calculations with the 
SLy4 force [23[ ■ The Coulomb potential Vc (r) is well de- 
fined physically and should not be folded. We have ap- 
proximated the Coulomb potential by Vc(r) = Z c v c (r), 
where v c (r) is the single-proton Coulomb potential given 
by the mean-field calculation. In the spherical case, the 
folding factor A is determined with the Bohr-Sommcrfcld 
quantization condition [Til ]. 
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^\Q* L -V(r)\ = (2N+l)^ = (G-L+l)^ (3) 



where r\, ri (and r% in Eq.(5)) are the classical turning 
points obtained by V(r) = Q* L (the decay energy). The 
global quantum number G is estimated by the Wildcr- 
muth rule, depending on the configurations of valance 
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TABLE I: The calculated a-decay widths of excited states 
in 8 Be and 20 Ne, with comparison with the experimental 
data [H. 





E*, 

J 


(G,L) 


L Q (calc.) 


r a (expt.) 




[MeV] 




[keVl 


fkeVl 

L J 






8 Be (Qn= 


=91.8 keV) 




0+ 


0.00 


f4 0) 


7.8 ev 


6.8±1.7 ev 


2+ 


3.04 


(4, 2) 


1.6 xlO 3 


1.5xl0 3 


4+ 


11.40 


(4, 4) 


2.3xl0 3 


?«3.5xl0 3 






20 Ne (Q = 


-4.729 MeV) 




6+ 


8.78 


(8, 6) 


0.17 


0.11±0.02 


8+ 


11.95 


(8, 8) 


90 ev 


35±10 ev 


1" 


5.79 


(9, 1) 


18 ev 


28±3 ev 


3" 


7.16 


(9, 3) 


4.7 


8.2±0.3 


5" 


10.26 


(9, 5) 


54 


145±40 


~r 


15.37 


(9, 7) 


120 


110±10 


9" 


22.87 


(9, 9) 


116 


225±40 



nucleons [24| . The cluster decay energy from an excited 
state is given by, 

Q* L = Qo + E*, (4) 

where Qq is the decay energy from the ground state, and 
Ej is the excitation energy of a given state with the spin 
J. The decay process can occur only if the state has a 
positive Q* L value. 

The partial decay widths can be calculated by [l4| , 

drk(r)} 

1=^ jn , (5) 

/ dr/2k(r) 

J ri 



where k(r) = y j^\Q*l ~ V( r )\ 1S the wave number, and 
P is the preformation factor of the cluster. For even- 
even nuclei, the extreme P=l assumption under using 
the Bohr-Sommerfeld condition can well reproduce the 
experimental half-lives of various cluster decays [TJ, H3] ■ 
The decay half- life can be obtained by T^/ 2 = H\n2/T. 

In the axially deformed case, the decay width can be 
approximated by averaging widths at various directions 
of the space as follows [25| 

,7r/2 

r = / r(0) sin(0)d0. (6) 

Jo 

To calculate the width T(9), a deformed cluster po- 
tential has to be employed. In the present work, the 
deformed potential V(r, 9) is constructed with axially 
deformed single-particle potentials w„(r, z), v p (r, z) and 
v c (r, z) that are given by the shape-constrained Skyrme- 
Hartree-Fock calculation [26| , with the folding factor de- 
termined at the spherical case. The variables (r, z) are 



TABLE II: The calculated a-decay widths for the = 0^ 
cluster band in 10 Be (Qo — — 7.413 MeV) at the different cases 
of the spherical and deformed shapes. 



r 


E}[MeV] 


r Q (sph.) 


r Q (def.) 


r a (z-axis) 


T a (expt.) 


o 2 + 


6.18 










2+ 


7.54 


0.51 ev 


0.85 ev 


4.1 ev 


22±7 ev [4] 


4+ 


10.15 


35 keV 


97 keV 


584 keV 


130 keV [5] 



TABLE III: The calculated 6 He-widths of the states belong- 
ing to the K n = 0^ band in 12 Be (Q o =-10.11 MeV). The 
AMD predictions are also given for comparison. 



r 


£}[MeV] 


r6 Hc (present) 


r 6Hc [9j 






[keV] 


[keV] 


o 3 + 


10.9 


410 


700 


2+ 


11.3 


285 


1 


4+ 


13.2 


190 


7 


6+ 


16.1 


34 


16 


8+ 


20.9 


3.7 





the cylindrical coordinates, and 9 is the angle between 
the symmetry axis and the radius for the cluster emis- 
sion. 

The most well-established cluster structures in light 
nuclei are in 8 Be and 20 Ne (24|, with an a-particle cou- 
pling to the magic cores of 4 He and 16 O, respectively. 
Their cluster structures can be well described by the 
semiclassical cluster model [H,[27]j. For 8 Be, we take G=4 
for the ground-state band, according to the Wildermuth 
condition For 20 Ne, we take G=8 for the K n =0 + 

band (the ground-state band) and G=9 for the K w =0~ 
band, as discussed in Refs. [UH^I- The calculated re- 
sults are given in Table I. It can be seen that the present 
calculations agree well with the experimental widths [28[ 
within a factor of 3. 20 Nc is a particularly interesting 
nucleus that can have very different structures for differ- 
ent states. While the = 0~ band with the sequence 
of = 1 _ ,3 _ ,...,9 _ has an almost pure a+ le O clus- 
ter structure, the + ground-state band has a consider- 
able mixture of the cluster structure and the deformed 
mean- field structure [2l|. For the 0~ band, both experi- 
ments and our calculations give remarkable large a-decay 
widths comparable with the Wigner limit [2^ | , indicating 
the significant a+ le O cluster structure. 

For heavier beryllium isotopes, the 2a-cluster struc- 
tures play an important role. For example, the well- 
known parity inversion in 11 Be is related to a large pro- 
late deformation of 10 Be due to a developed cluster struc- 
ture 0. The deformed mean- field calculation for 11 Be 
with the mixture of a sd neutron orbit is insufficient to 
gain enough binding for the J 71 = l/2 + state (26|. De- 
formations in beryllium isotopes arc dependent on the 
distance of the two a particles. In 10 Be, the cluster band 
with 0^(6.18 MeV), 2+(7.54 MeV) and 4+ (10.15 MeV) 
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members have been established experimentally 0, [3(| ■ It 
was found that a channel radius as large as 8±1 fm has to 
be adopted to reproduce the a-decay width for the 7.54 
MeV state [H, while the a- width of the 10.15 MeV state 
is associated with a smaller channel radius of 5^6 fm p|. 

From the shell-model viewpoint, the lo Bc(0f) struc- 
ture has a (sd) 2 configuration due to the large prolate 
deformation [3 11 ] . Hence we assume a global number G=6 
for the decay calculation. The calculated widths in 10 Be 
are displayed in Table II. It can be seen that calculations 
in the spherical case give much smaller a-decay widths 
compared with experiments. As predicted by the molec- 
ular orbital model [3l|, the a-a distance in the second 0+ 
state is about 4 fm. This corresponds to an axis ratio of 
2.5:1 (or a prolate deformation of /?2 ~ 1-1). With such 
a large deformation employed in Eq.(6), the a- width for 
the 10.15 MeV state can be reproduced reasonably. We 
found that the a-width of the 2+ (7.54 MeV) state is sen- 
sitive to the decay energy rather than the deformation of 
the cluster potential (the width can be reproduced with 
an increase of only 50 keV in the decay energy). This 
can be understood considering that, for a state near the 
threshold, the cluster tunnelling with a very small decay 
energy is dominated by the long tail of Coulomb potential 
(the tail is not sensitive to the deformation). 

Considering the a:2n:a molecular structure of the Of 
band we estimated the a- width by assuming that 
the a particle emits along the z axis with the fii = 1.1 
deformed potential. The calculated widths are given in 
Table II. It can be seen that the a- width for the 7.54 
MeV state is improved significantly. However, the a- 
width for the 4 + (10.15 MeV) state is overestimated. It 
was pointed out that the Of band states could contain 
the significant mixture of the 5 He+ 5 He configuration (3ll . 
[32j| . The mixture is particularly remarkable for the 10.15 
MeV state because it is almost on the 5 He+ 5 He decay 
threshold, which could result in a reduced a-decay width. 

For the neutron-rich 12 Bc nucleus, some resonance 
states were recently observed, decaying to a+ 8 He and 
6 He+ 6 He channels [1, 0]- A rotational band with a 
6 He+ 6 He structure has been suggested to have a very 
large momentum of inertia (h 2 /2^ = 0.14 MeV) Q. The 
AMD calculations interpreted that resonances happen in 
molecular states built on the O3 configuration Con- 
sidering that the highest state of this band would be the 
J 11 = 8+ member % Q, we take G=8 in the 6 He+ 6 He 
decay calculation. A spherical cluster potential is as- 
sumed. In Table III, the obtained widths are compared 
with the AMD predications [§]. As stated in [J, the 
6 He-width of the Of state is very large due to the lack of 
centrifugal barrier. It is shown that the present widths 
significantly decrease with increasing angular momenta. 
In these states, a+ 8 He decays have also been observed Q 
but their description is beyond the cluster tunnelling 
model because of too large decay energies. 

Carbon isotopes heavier than 12 C can have three- 
centre cluster structures with different geometric shapes. 
As suggested by von Oertzen [13, HH, there are possi- 
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Prolate Oblate 

FIG. 1: The schematic picture for the two modes of a decays 
in carbon isotopes, corresponding to the prolate(linear) and 
oblate(triangle) structures, respectively. 



ble prolate(lincar) and oblate (triangular) cluster states 
in carbon isotopes. The K* = 0+ and 3j~ bands in 12 C 
are considered to have triangular configurations. The 
0^ state of 12 C has been predicted to be an analogue 
of Bose condensate [Hj]. Recently, a decays from the 
excited states of 13 C and 14 C have attracted intense ex- 
perimental studies (see e.g., [H, EH)- 

The oblate and prolate configurations certainly lead 
to different calculated decay widths. As shown in Fig.l, 
the cluster potential radius in the oblate case(# = 7r/2) 
is smaller than in the prolate case(6> = 0). Therefore, 
decays in oblate case would be suppressed. In 12 C, the 
4+ (14.08 MeV) and 3^(9.64 MeV) states belong to the 
= Of and 3j~ oblate bands, respectively. The de- 
cay channels of these low-lying states arc nearly pure 
a decays(r Q ~ T) [H, HH, providing a good chance to 
study the deformation effect. In Table IV, with a de- 
formed cluster potential of /?2 = 0.8 (corresponding to 
an axis ratio of 2:1), the obtained a-widths of the 3~ 
and 4 + states at 6 = tt/2 direction (the oblate case) 
agree well with the experimental widths. However, their 
a-widths at 9 = direction (the linear case) are calcu- 
lated to be 243 keV for the 3" state and 2.2 MeV for 
the 4 + state, which are about one order of magnitude 
larger than the oblate calculations. This indicates that a 
decays are strongly hindered in triangular configurations 
compared with linear configurations. 

The a decay from the Of state in 12 C has been studied 
with the triple cluster model [34] . This state has a decay 
energy of 288 keV to 8 Be(g.s.)+a and an energy of 380 
keV to 3a particles. For the decay to 8 Be+a, the present 
calculation with a linear structure (P2 = 0.8) gives an a- 
width of 5.9 ev that agree with the observed value of 8.5 
cv. The calculations assuming a triangular configuration 
leads to a small width of 1.6 ev. For the decay to the 3a 
particles with a triangular structure, our model estimates 
a width of 36 ev that is close to the predictions of the 
AMD [35| and the triple cluster model [34 [. However, the 
contribution of the direct 3a-decay process to the total 
width is less than 4% obtained in the experiment [3?| . 
Hence, the linear cluster decay to 8 Be+a should domi- 
nate the decay of the Of state in 12 C. 
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TABLE IV: The calculated a-widths of excited states in car- 
bon isotopes. The experimental total decay width T are taken 
from HH for 12 ' 13 C and from [HI for 14 C. 



J" 


E}[MeV 


(G,L) 


L a (calc.) 


T(expt.) 








[keV] 


[keV] 




12 C (Q =-7.366 MeV) 




■ 1 


7.654 


(6, 0) 


5.9 ev 


8.5 ev 


37 

1 


9.641 


(5, 3) 


17 


34 


4+ 

1 


14.08 


(4, 4) 


158 


258 




13 C (Q = -10.647 MeV) 




K n = 


3/2" 








5/2" 


10.82 


(6, 2) 


0.2 ev 


24 


7/2" 


12.44 


(6, 2) 


84 


140 


9/2~ 


14.13 


(6, 4) 


32 


150 


11/2 


16.08 


(6, 4) 


316 


150 


if* = 


3/2+ 








3/2+ 


11.08 


(7, 1) 


41 ev 


<4 


5/2+ 


11.95 


(7, 1) 


100 


500 


7/2+ 


13.41 


(7, 3) 


112 


35 


9/2+ 


15.28 


(7, 3) 


888 




11/2 


+ 16.95 


(7, 5) 


125 


330 




14 C (Q = -12.011 MeV) 




5" 


14.87 


(7, 5) 


0.4 


35 


6+ 


16.43 


(6, 6) 


0.12 


35 


;r 


12.58 


(7, 3) 


22 ev 


95 


5" 


15.18 


(7, 5) 


31 


50 


7" 


18.03 


(7, 7) 


14 


70 


6+ 


14.67 


(6, 6) 


0.21 


40 


TABLE V: The calculated a-widths of the 14 C excited states 


decaying to 10 Be(g.s.^ 


+a and 10 Be(2 


+ )+a channels. 




18.5 MeV 


19.8 MeV 


r 


10 Be(g.s.) 


10 Be(2+) 


10 Be(g.s.) 


10 Be(2+) 


4+ 


1.2 MeV 


0.8 MeV 


2.7 MeV 




5" 


0.28 MeV 


0.39 MeV 


0.7 MeV 


1.2 MeV 


6+ 


11 keV 


30 keV 


41 keV 


230 keV 


7" 


0.8 keV 


2.3 keV 


3.5 keV 


25 keV 



For 13 C, the K* = 3/2 ± bands constructed on the 
a+ 9 Be[3/2", g.s.) structure were suggested experimen- 
tally [Ty] ■ The parity doublet bands are related to the 
reflection asymmetric chain configurations, correspond- 
ing to a very large prolate deformation (h 2 /2^s = 0.19 
MeV) [1(3]. We assume that the core (i.e., 9 Be) has a 
similar deformation to that of 10 Be with @2 ~ 0.6 exper- 
imentally [37[- In the a- width calculations, we take G=6 
and 7 for the = 3/2" and 3/2+ bands, respectively, 
considering the different parity of the bands. The angular 
momenta L that the a particle carries are obtained by the 
angular momentum selection among the parent, daughter 



and a particle. Furthermore, for the odd nucleus 13 C, we 
take the a preformation factor P=0.6 that was adopted in 
the systematical a-decay calculations of odd nuclei [HI]. 
The calculated a-widths of the K v = 3/2+ 1 band states 
in 13 C are listed in Table IV. The experimental a- widths 
have not been available, but we give the experimental 
total widths T [l^, [28[ of the states for comparison. 

The calculated a-widths for 13 C are in general smaller 
than experimental total widths as expected. How- 
ever, calculated values for the 11/2~(16.08 MeV) and 
7/2+(13.41 MeV) states I1JJ are larger than experimen- 
tal total widths. In Ref. [28| , the two states were assigned 
with (7/2+) and (9/2~), respectively. The recent experi- 
ment [l2| suggests that the 13.41 MeV state is connected 
to an oblate structure and the 16.08 MeV state has a 
positive parity. Calculations can be changed significantly 
due to the different assignments of spins. 

The predicted a-widths for 14 C decaying to the ground 
state of 10 Be are also presented in Table IV. The bands 
built on the 0^(6.589 MeV) and 3^(9.801 MeV) states 
were suggested to have oblate (triangular) structures, 
while bands built on the 0^(9.746 MeV) and 1"(11.395 
MeV) states have prolate (chain) structures [ll[. In Ta- 
ble IV, the 5" (14.87 MeV) and 6+ (16.43 MeV) states 
have oblate structures and other states are prolate [Til ] . 
The calculated width of the 3" (12.58 MeV) state is 
remarkably smaller than the experimental total width. 
This would indicate that the neutron emission chan- 
nel dominates the decay of the state that is near the 
a-decay threshold. The decay channel of 14 C(16.43 
MeV)^ 10 Be(2+,3.37MeV)+a has also been observed 
though it is very weak [H[ . The decay width is estimated 
to be 5 ev, giving a small branching ratio compared with 
the 10 Bc(g.s.)-|-a decay channel. 

It is interesting that, in 14 C, the 18.5 MeV and 19.8 
MeV states have strong channels decaying to both the 
ground state and the first 2+ state of 10 Be Q3, [Hj]. 
No experimental assignments of spins and parties have 
been available for the two states of 14 C. We estimated 
their a-widths with assuming J 71 " = 5~,7~(G=7) and 
J n = 4+,6+(G=6). The calculations are preformed in 
the simple spherical shape. As shown in Table V, the 
calculated a-widths are sensitive to the assignments of 
angular momenta. With the spin of 4~5, broad widths 
are obtained for both channels as observed in experi- 
ments [Hill. 

In summary, we have calculated decay widths for the 
interesting cluster states in beryllium and carbon iso- 
topes, such as the a+ 6 He, 6 He+ 6 He and 8,9 ' 10 Be+a sys- 
tems, using a uniform folding cluster potential that has 
good descriptions for various cluster radioactivities in 
heavy nuclei. It has been shown that the model can well 
reproduce the a-widths of 8 Bc, 12 C and 20 Nc. The dis- 
crepant decay properties of the 7.54 MeV and 10.15 MeV 
states in 10 Be have been discussed. The 6 He-widths of 
possible molecular states in 12 Be and a-widths of excited 
states in 13 < 14 C are predicted. 
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